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Abstract 

The quantum dynamics of linear molecular aggregates in the presence of Sq— ^-Si and 
Sq— 5-82 transitions is investigated putting emphasis on the interplay between local non-adiabatic 
S2 to Si deactivation and Frenkel exciton transfer. The theoretical approach combines aspects 
of the linear vibronic coupling and Frenkel exciton models. Dynamics calculations aie per- 
formed for the absorption spectrum and the electronic state populations using the multiconfig- 
uration time-dependent Hartree approach. As an application perylene bisimde J-type dimer and 
trimer aggregates are considered, including four tuning and one coupling mode per monomer. 
This leads to a dynamical model comprising up to 7 electronic states and 15 vibrational modes. 
The unknown non-adiabatic coupling strength is treated as a parameter that is chosen in accor- 
dance with available absorption spectra. This leaves some flexibihty that can be limited by the 
clearly distinguishable population dynamics. 
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Introduction 



Exciton- vibrational coupling (EVC) plays an important role in the dynamics and spectroscopy of 
natural and artificial molecular aggregates.-"- In particular the non-perturbative and non-Markovian 
regime has attracted recent interest and a number of approaches have been developed to cope with 
such situations.-^— Often the signatures of strong EVC are visible as vibrational side bands already 
in linear absorption or emission spectra.—""— However, EVC is not the only manifestation of the 
interaction between electronic and nuclear degrees of freedom. A frequent theme in the dynamics 
of excited electronic states of molecules are non-adiabatic transitions signifying the breakdown of 
the Bom-Oppenheimer approximation.— The topic of this contribution is the interplay between 
EVC in exciton transfer and local non-adiabatic transitions. This issue has been intensively stud- 
ied in the context of exciton-exciton annihilation. There highly excited states at energies about 
twice the Sq— ^-Si excitation energy are involved. In contrast we focus on situations where the 
monomers comprising the aggregate have overlapping electronic transitions. In terms of exciton 
theory this will lead to the emergence of two excitonic bands being connected by interband transi- 
tions which are due to non-adiabatic couplings. This situation did not receive much attention so far, 
although it should be rather common for many dyes. For instance, Kobayashi and coworkers have 
made vibronic coupling (Herzberg-Teller coupling) between Q and B state derived exciton man- 
ifolds responsible for observed vibrational dynamics after Q-band excitation in porphyrin type J- 
aggregates.— Recently, Q- and B-band dynamics of bisporphyrin have been stduied using transient 
absorption and two-dimensional spectroscopy.— Transient absorption data revealed time scales of 
associated intraband relaxation and S2-S1 internal conversion in the few hundred femtosecond 
regime. Interband transitions have also been discussed for tubular aggregates. -'^^ However, in that 
case the excitonic manifolds correspond to different structural elements such as inner and outer 
tubes which have some residual electronic coupling. As a note of caution we stress that the present 
non-adiabatic coupling between different local electronic states should not be confused with the 
discussion of non-adiabaticity in the context of the potential energy surfaces (PES) coupled via the 
Coulomb interaction between different monomers.-^^s^ii^Si This also holds true for studies of molec- 
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Figure 1 : (color online) Schematic view of the topology of the diabatic PES for a molecular dimer 
(site index m = 1,2) including one vibrational coordinate and three electronic states, Sg-Sa, per 
monomer. The different excitation states of the dimer are indicated. 

ular complexes that are not treated as an aggregate but a supermolecule. For example, Leutwyler 
and coworkers interpreted splittings at the electronic band origin of Sq— >Si and Sq— 7'S2 transitions 
of hydrogen-bonded homo-dimers in terms of excitonic coupling quenched by nuclear tunneling.— 
The situation we will focus on is illustrated in Fig. [T]for the case of a dimer with one harmonic 
vibrational coordinate and three electronic states for each monomer. The lowest state, (So, So), 
corresponds to the situation where both monomers are in their electronic ground state, Sq. Next 
in energy there are two diabatic PES corresponding to the situation where a single Sq-Si type 
excitation is present, either in monomer 1 (Si, So) or 2 (So,Si). For the depicted symmetric case 
the two PES intersect along the diagonal of the plot. In the adiabatic representation one has an 
avoided crossing due to the Coulomb coupling (cf. ref. ll^. If there is a close lying higher in- 
tramonomer electronic state, S2, another set of Coulomb coupled diabatic PES will be present 
corresponding to the excitation states (S2,So) and (So,S2). Since the local shifts of the electron- 
ically excited states along the considered nuclear coordinate in general are different, the related 
PES will intersect. Assuming that besides the depicted tuning modes, there are some coupling 
modes, conical intersections (CIs) will emerge along the PES intersection lines.— The two types 
of state couplings will manifest themselves in the absorption spectral lines shapes. Furthermore, 
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one might think of interesting quantum dynamical implications. For instance, consider the limiting 
case where a system initially prepared in a local S2 excitation will evolve either in the S2 manifold 
or after non-adiabatic transition in the S 1 manifold. Which dynamics is relevant will depend on the 
ratio of non-adiabatic and Coulomb coupling. In fact if both are of the same magnitude intra- and 
interband transfer may be strongly mixed. An important aspect for the dynamics is the coupling 
between electronic and nuclear degrees of freedom because it influences both local non-adiabatic 
transitions and transfer (e.g. shifts of PES). In general, it will not be a single vibrational mode that 
mediates both processes. 

In the present contribution the scenario of Fig. [T]will be studied for the exemplary case of a 
modified perylene bisimide (PBI) dye forming J-aggregates in solution.— The choice of the system 
is motivated by our previous studies, which focussed on the quantum chemical and spectroscopic 
characterization of monomers and small aggregates.-2^>22 Here we will extend our model to include 
the S2 state as well as non-adiabatic 82—7-81 transitions. In the following section we start with 
outlining the theoretical approach which combines the Frenkel exciton Hamiltonian with vibronic 
coupling theory in the diabatic representation. Next some details of the quantum dynamics method 
for calculating absorption spectra and population dynamics are given. This section concludes with 
a summary of the model parameters. The Results section focusses on the spectral features due to 
the different couplings as well as on the corresponding regimes of population dynamics. 

Theoretical Methods 

Frenkel-Exciton Vibronic- Coupling Hamiltonian 

In the following we will give a short summary of the Frenkel exciton Hamiltonian with special 
emphasis on the vibronic coupling between different local diabatic electronic states, \(Pm,a) (m = 
1 . . .A/^agg is the site index and a = 0, 1,2 counts the electronic states). The aggregate Hamiltonian 
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can be written as— 



)\(pmM(pn,b){(Pn,c(Pm,d\- (1) 

m ab innabcd 

Within the shifted harmonic oscillator model the diagonal on-site part is given by electronic 
energy) 

\ TJ , Y^ fi(Om.a{^) ( , ^2 \ 

Um,a ~ E,na + £^Wm,a(<^)^,n,a(^)e„,,^. (3) 

Here {2m ,e} denotes the set of dimensionless electronic ground state normal modes at site m, 
(X>m.a{^) are the respective frequencies, and gm.a{^) is the linear coupling in state |(Pm,a) (8m,a{^) = 
{dUm,a/dQ,n,^)/ho),n.a{^))- The coupliug of a particular mode to the electronic transition can be 
characterized by the (dimensionless) Huang-Rhys factor Sm,a{^) = ,?m,a(^)/2- These modes are 
also called tuning modes since they influence the energy gap between different electronic states. 

For the off-diagonal elements a linear expansion in terms of the set of coupling modes {Qm.ri} 
is assumed to be valid (linear vibronic coupling (LVC) model), i.e.— 

Hm{a,b) =Y,^m,-nia,b)Q„i,T^ (4) 

with X,n,r\{a,b) being the respective coupling matrix. The A^vib modes of the mth site will be 
comprised into the vector Q„, = {Qm,^}) 

Electronic transitions at different sites are coupled via the Coulomb matrix element Jmn{cimb,n Cndm) 
in Eq. ([T])- In the following we will assume that concerning the Coulomb coupling only close to 
resonant interactions involving ground state transitions need to be accounted for, i.e. the only non- 
zero elements of Jmni^mbm^-ndni) are Jmni^ m^m ^n^ni) JmniS^m^m^n^m) <ind Juifii^irfim'^n^m) 

Jmn{0,n'2ni 0„2,„) . Further, J^n will be assumed to be independent of the nuclear coordinates. 

Electronic excitations are usually classified according to zero-, one-, two- etc. excitation states. 
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Restricting to the one-exciton space we have the following completeness relation 



l=|0)(0|+£ £ \ma){ma\, (5) 

m a=l,2 

with 

io) = n I *P"^'0) ' i'"^) = I'P"'.") n I ^kfi) ■ (6) 

Thus we have the correspondences: |0) — (So, So, ...), I la) — ^ (Sq, Sq, ...) and so on. The 
aggregate Hamiltonian expressed in this basis can be written as 

//agg = (7) 

= £//„(0,0)|0)(0|=^o|0)(0|, (8) 

m 

//(I) = £ £ [5„„(4z,(^o + t/,.,a) + (l-5a^)//™(a,^)) 

mna,fo=l,2 

+ 7m)i('3m,0„,Z7„,0m)]|ma)(n^| . (9) 



Quantum Dynamics 

The time-dependent Schrodinger equation will be solved employing the multiconfiguration time- 
dependent Hartree (MCTDH) method.-^^^ To this end the state vector is expanded in terms of the 
diabatic basis, i.e 

imO) = Eza(Q;0|a), (10) 

a 

where a runs over all electronic configurations (cf. Eq. (??)) and the nuclear coordinates are 
comprised in the D = N^gg x A'vib dimensional vector Q = (Qi, . . . , Qw^gg) = • • • , 2d)- 
The nuclear wave function for each diabatic state is expanded into MCTDH form as follows 

"jV-"jD 

Jv-Jd 
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Here, the Cj^ j^{t) are the time-dependent expansion coefficients weighting the contributions of 

(a) 

the different Hartree products, which are composed of single particle functions (SPFs), [Qk'j), 
for the kih. degree of freedom in state a. 

The aggregate models will be characterized by means of the absorption spectrum. Here, a time- 
dependent formulation will be used, i.e. the absorption spectrum at T = K is expressed as^ (/q 
normalization constant) 

Iico)=Io(0 E Re / dte''''-^'/''^'{'¥o\S'^e-'"^"^^'/^S'^\'¥Q), (12) 
where d^'^ is the z-th directional component of the dipole operator. For the latter we assume the 



Condon approximation with constant S 
be valid. Notice that according to ref. 
polarization, i.e. 



— )-Si and Sq— )'S2 transition dipole moments d,„^a=i,2 to 
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the Sq— J-Si and the Sq— )'S2 transitions have different 



= £(J„,,i|mi)(0|+h.c.), (13) 
111 

d^y^ = £(J,„,2|m2)(0|+h.c.). (14) 

m 

Further, in Eq. (??) T,- is a parameter mimicking the finite line width (dephasing time) of the real 
system. Note that we will use different values for the Sq— J-Si and 80^82 transitions in order 
to obtain better agreement with experiment. This is motivated by the fact that in contrast to the 
orbitals involved in the Sq— !>Si transition those of the 80—^82 transitions are partly localized on 



the phenyl bay substituents (cf. ref. 



261) . Their Soppiness will cause stronger fluctuations of the 



transition energy in a solvent environment. 



Computational Methods 

The setup for the electronic structure calculations has been detailed in refs. 
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and 



27 



In brief 



the gas phase ground state equilibrium geometry of the PBI monomer has been obtained using 
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density functional theory (DFT) with the B3LYP functional and a 6-31 IG* split- valence basis set. 
The first two electronically excited singlet states were calculated employing time-dependent DFT 
(TDDFT). Harmonic analysis of the ground state vibrations and projection of the forces at the 
vertical transitions have been used to obtain frequencies and Huang-Rhys factors. All quantum 
chemistry calculations have been performed with the TURBOMOLE program package.— For the 
Coulomb coupling between the Sq— t-S i transitions we have used a value which is in accord with the 
previous calculation— and an estimate from experiment (see below).— For the Sq— )-S2 transitions 
we have scaled the Sq— J-Si value according to the different calculated transition dipole moments 
assuming the dipole approximation for simplicity. 

Wave packet propagations have been performed using the Heidelberg program package.— 
Three systems have been considered, i.e. monomer, dimer, and trimer having 3, 5, and 7 dia- 
batic states, respectively. For each monomer five nuclear degrees of freedom are included into 
the model with parameters to be specified below. Four of these, Qm.i-Qm,4, act as tuning modes 
for the system. The fifth mode, Qm.c, represents a coupling mode. As primitive basis we have 
chosen a 35 point harmonic oscillator discrete variable representation in the interval [—7.5 : 7.5] 
and [—6.5 : 6.5] for the low- and high-frequency tuning mode, respectively. For the coupling mode 
we have used 35 points in the interval [—6 : 6]. For the single SPF basis the multi set method was 
used (cf. Eq. (??)). Further, we employed mode combination with respect to the tuning modes, i.e. 
alike modes of the different monomers (m = 1, . . . ,A^agg; ^ = 1, •••4) have been combined to 
four A^agg-dimensional sets of basis functions. As convergence criteria for the SPF basis we used 
that the maximal population of the least occupied SPF should be below 10 in the monomer and 
dimer case and below 5 ■ 10 in the trimer case. This required from 9/2 to 25/6 SPFs per state and 
combined tuning/coupling modes for the weak and strong coupling case, respectively. 

Absorption spectra have been calculated via two independent wave packet propagations ac- 
cording to the sum in Eq. (??) using the transition specific dipole operator in Eq. (fT3T ). Diabatic 
state population dynamics is investigated starting from two different initial conditions: (i) exci- 
tation from the So to all S2 states according to the dipole operator d^^^ and (ii) excitation of the 
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first monomer m = 1 only. The latter condition, although being artificial in terms of experimental 
realization, is better suited for investigation of the interplay between local non-adiabatic transitions 
and energy transfer. 



Results 



Application to PBI 

The electronic parameters have been adopted such as to present the situation in the J-aggregate 



forming PBI dye (see, refs. 
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- 27|) . Thereby we have neglected a possible heterogeneity and 
assumed that all monomers have identical properties and Coulomb couplings between neigh- 
boring monomers are the same throughout the aggregate. Specifically, we have chosen E,n,i = 
El = 2.13 eV and 2 = £2 = 2.1 A eV. For the Coulomb couplings we use the following values: 
Jmni^mOn, ^nOm) = —500 cm^^ and 7,„„(2„,0„, 2„0m) = —150 cm^^ Since we will show normal- 
ized spectra, only the ratio between the transition dipole moments is of relevance. According to 



ref, 
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this ratio is given by d,n^\/d,n,2 = 1-85. 
In Fig. [2] we present the Huang-Rhys factors for the monomeric Sq— t-Si and Sq— 7'S2 transitions. 
First we notice that the contributions in the low-frequency region (below 100 cm^) are of rather 
high magnitude. Given the potential problems of this spectral region, e.g. with respect to necessary 
anharmonic corrections, we will not consider frequencies below 100 cm^^ In order to set up a 
simple model, which is extendable to larger aggregates, it will be necessary to comprise the many 
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we could 



modes with appreciable Huang-Rhys factors into some effective modes .-^^^^ In ref. 
obtain a reasonable fit to the linear Sq— !>Si absorption spectrum by combining all modes above 
800 cm^^ into one effective mode. Here we will introduce a second effective mode for the range 
between 100 and 800 cm ^ The resulting effective mode parameters for both transitions are 
compiled in Tab. [T] From Tab. [T]we notice that in accordance with Fig. [2] the two excited state 
configurations have nuclear displacements in different directions. 

In the following we will consider three different models, i.e. the cases Nagg = 1, 2, 3. In each 




1000 1500 
wavenumber [cm-i] 



2000 



Figure 2: Huang-Rhys factors for the Sq— !>Si and 80^82 electronic transitions of a PBI monomer. 
Notice that the low-frequency part has been cut-off since for modes in this region the harmonic 
approximation is less reliable. 



case the monomer is described by three electronic states, four tuning, and one coupling mode. 
Conceming the tuning modes we have used the effective mode model of Tab. \T\ Since there is no 
information with respect to the coupling mode we have assumed that its frequency is 1000 cm^^ 
and the LVC constant A is treated as an adjustable parameter. Note that the qualitative behaviour 
is not much influenced by the actual choice of the frequency as has been confirmed for a value of 
100 cm^^ . In order to have some guidance for the choice of A we will compare linear absorption 
spectra with experimental results.—"^ Since this allows for some ambiguity we will further inspect 
the population dynamics of diabatic states, which in principle could be obtained, e.g., by fitting 
from pump-probe spectra.- 

In order to gain some insight into the topology of the PES we have plotted in Fig. [3] different 
cuts for a dimer in adiabatic representation and using the antisymmetric combination of monomer 
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Figure 3: (color online) (a) Adiabatic PES of the Si (blue) and S2 (red) derived states for the dimer 
along antisymmetric combined modes = {Qi_^ — Qi,^)/^- CI's exist between the lower 
state of the S2 band and the higher state of the S 1 band for non- vanishing coupling matrix element 
within the LVC model. Panels (b) and (c) show a cut along the Q2 and Q4 mode, respectively, 
for better visibility of the energetics. 



coordinates, = —Qi.^)/"^^- Panel (a) shows a global view whereas panel (b) and 
(c) present a cut along that high-frequency mode which is displaced in the Si and S2 state, respec- 
tively. Important for the subsequent discussion are the Coulomb coupling induced splittings around 
= and the region of approach of Si and S2 surfaces for | > 4. Here, CIs emerge due to 
the LVC when moving along the orthogonal mode Qm,c- The energetic position of the CI depends 
mainly on the difference of Huang-Rhys-factors for the Si and S2 states, but is also influenced 
by the band splitting due to the Coulomb coupling. It should be noted, that the initial dynamics 
after photo-excitation proceeds not along this coordinate, but along the symmetric combination 
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Table 1 : Effective mode parameters which are calculated from the full set of modes (denoted by a 
tilde) as follows: frequency (Om^ai^) = l/S„i,ai^)'£^'S,n^aiQa},niQ (given in cm^^), Huang-Rhys 
factor S,„^ai^) = L^''^m.a(C)' whcrc the prime should remind on the definition of the frequency 
intervals for the two effective modes (see text). Note that no excited state frequencies are calculated 
in the original full set, i.e. only information about the vertical transition geometry is included. 



mode (^) 








1 


297 


0.71 





2 


1416 


0.44 





3 


400 





0.49 


4 


1354 





0.47 



= (2i,<^ + 22,1^ ) / v^- Hence, for the dimer reaching the CI requires a backscattering from the 
repulsive potential hit upon moving along Q^^ (cf. Fig.[T]). 



Absorption Spectra 

In Fig. |4]we compare the absorption spectra of the different systems for two values of the LVC 
strength A with experimental data obtained for monomeric— and aggregated-^- FBI. The agree- 
ment is fairly good, given the approximate nature of the present reduced dimensionality model as 
well as possible errors due to the quantum chemical method. Further note that in experiment the 
aggregates are most likely longer than the trimer considered here, although it has been argued that 
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the effective coherence length of the exciton is two monomers only; see also discussion in ref 
The effect of J-aggregation on the spectrum is seen, independent of the coupling strength A, i.e. the 
red-most band of the spectrum increases in intensity relative to the higher energetic part. Because 
of the higher Coulomb coupling strength between the S i states as compared with the S2 states the 
red shift is stronger in the So^>Si band than in the Sq— )'S2 band of the spectra. As expected this ef- 
fect of concentration of oscillator strength is more pronounced for the trimer as compared with the 
dimer. Although the experimental spectrum is considerably broadened it is instructive to inspect 
a high-resolution simulation, which gives insight into the underlying transitions. For the case of 
the monomer (panels (a) and (b)) the Sq— J-Si and Sq— )'S2 derived Franck-Condon progression are 
clearly discernible. They partly overlap in the range of the Sq— ;'S2 0-0 transition around 2.7 eV. The 
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Figure 4: (color online) Absorption spectra of monomer (a,b), dimer (c,d), and trimer (e,f) for 
two different values of the LVC constant A. The comparison with the experimental data (ref. 
26l]35L green line) is made using two different resolutions (blue line = 30 fs (Sq— J-Si ) and 11 fs 
(So^S2 ), red line Ty = 400 fs). 



effect of increasing the LVC constant is most notable in the Sq— )'S2 band although some oscillator 
strength redistribution is observed even for the low-energy part (2.1-2.3 eV) of the Sq— j-Si band. 
In case of the dimer (Fig. |4tc,d)) the Sq— J-Si part of the spectrum is less affected by the LVC and 
in fact the effect of increasing the aggregate to a trimer (Fig. |4te,f)) has a more pronounced in- 
fluence. Clearly the Sq— J-Si derived band is shaped by the Coulomb coupling. In contrast for the 
Sq— band both LVC and Coulomb couplings play a role. Besides the mere magnitude of the 
couplings it is the energetic position of the CI which is responsible for this behavior (cf. Fig. |3]). 
Most important in this respect are those high-frequency modes with a large Huang-Rhys-factor for 



13 




50 100 150 200 250 300 350 400 50 100 160 200 250 300 350 400 

t lis] I lis] 



Figure 5: (color online) Diabatic population dynamics of S2 (red) and S 1 (blue) states for monomer 
(a,b), dimer (c,d), and trimer (e,f) and for two different values of the LVC constant X. The initial 
excitation is to the collective S2 band according to the dipole operator, Eq. ([T3T l. For the aggregate 
case the populations of the different locally excited diabatic states have been summed. Notice that 
due to symmetry all alike local excitation states behave in the same way. 

the S2 states, which give rise to low lying CI's. 

Summarizing Fig.|4]it can be stated that (i) given the approximations involved the experimental 
monomer and aggregate spectra are reasonably reproduced and (ii) based on these spectra alone, 
the magnitude of the LVC strength cannot be unambiguously determined. We note, however, that 
value beyond the used 0.1 eV lead to markedly different absorption spectra and therefore can be 
excluded. 
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Population Dynamics 

Whereas the low-resolution spectra do not show a pronounced dependence on the LVC strength 
over a reasonable parameter range, the population dynamics is, of course, rather sensitive to this 
parameter. This can be seen by comparing the left and right columns of Fig. [5l For the larger value 
A = 0. 1 eV the S2 to Si interband population transfer is almost complete within 400 fs for all cases. 
We further notice that oscillations are more pronounced for A = 0. 1 eV whereas for A = 0.025 eV 
the population decay appears to be almost exponential. Comparing the monomer with the two 
aggregates we observe a non-monotonic dependence on aggregate size for the smaller LVC. Here 
the interband decay becomes faster when going from the monomer to the dimer, but slows down if 
the aggregate is increased to a trimer. In contrast to the weaker LVC the final population after 400 
fs decreases with aggregate size for the stronger LVC. 

The observed behavior can be rationalized as being the net outcome of two opposite effects. 
First, the energetic position of the CI decreases with increasing aggregate size due to the influence 
of the Coulomb coupling (see Fig. [3]). Since this leads to a better resonance at the vertical transi- 
tion energy of the wave packet, the population transfer becomes faster as seen in the dimer case. 
However, as noted before the wave packet reaches the CI only after being backscattered from the 
repulsive wall in Qrefxi+ direction. This is in contrast to the monomer case where already the initial 
motion of the wave packet would be towards the CI. The broadening of the wave packet associ- 
ated with the backscattering causes a less effective interband transfer in the weak LVC case. The 
effect of backscattering is quantified in Fig. |6^ where we plot the time-averaged reduced density 
differences 

muQj) = [dt J dQ' (|^wj(Q,Ol'- l^woj(Q,Ol') , (15) 

where Qi and Qj denote selected coordinates, dQ' includes all other coordinates and the index 
"wJ" and "woJ" indicates a propagation with and without Coulomb interaction, respectively, up to 
final time T. 

Clearly, the presence of the Coulomb coupling leads on average to less density in the range 
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Figure 6: Wave packet dynamics in the S2 state of monomer 2 of the dimer model according to 
Fig. [5] averaged over a time interval of 400 fs. (a) Difference of diabatic densities Eq. (??) between 
simulation with (|^wj|^) and without (|^woj|^) Coulomb coupling (A = 0.025eV), (red: |^wj|^ > 
l^wojp. blue: l^wJp < l^wojp)- The horizontal line shows the crossing seam between S2 and Si 
states, (b) Difference of diabatic densities between simulation with small (|^2=o 025evl^) and large 

(|^A=o.ievl') LVC (red: |^A=o.ievl' > l^l=o.025evl^ blue: |^i=o.ievl' < l^l=0.025evp)- The 
black dot marks the CI geometry, (contours in interval [— 1 .2 : 1 .2] x 10 ^ (a) and [—7 : 1 .5] x 10 ^ 
(b) are chosen such as to promote visual clarity; wave functions are normalized to compensate for 
the population decay) 



of the crossing seam. In the dimer case the shift of the CI still dominates the effect of wave 
packet spreading, making the interband transfer faster as compared with the monomer. For the 
trimer however, the increased number of vibrational degrees of freedom allows for a pronounced 
intraband redistribution of the wave packet on the S2 manifold what slows down the interband 
decay. 

Finally, we comment on the observation of oscillations in the strong LVC case in Fig. |51 In 
case of weak coupling there is some continuous but relatively small leakage of density out of the 
S2 states, whenever the wave packet hits the crossing region. In contrast for strong LVC the S2 
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Figure 7: (color online) Population dynamics of local single excitation states (see key) for dimer 
(a,b) and trimer (c,d) and for two different values of the LVC constant X. The initial excitation is 
to the local S2 state of the first monomer {m= 1). 



population drops by about 40 % already on the first attempt. The subsequent oscillations are rather 
a signature of quantum interferences of parts the wave packet running on different states and being 
projected onto the diabatic electronic basis. The situation is visualized in Fig. [6b where we show 
the time-averaged density differences along two selected modes for the cases of small and large 
LVC in the dimer model. The wider spread of density for the larger LVC is clearly discernible. 
Finally, it should be noticed that these oscillations are likely to disappear upon increasing the di- 
mensionality of the model to approach a quasi-continuous spectrum of exciton- vibrational states. 

Deactivation versus Transfer 

In the following we will investigate the interplay between transfer and local deactivation in more 
detail. To this end we have chosen an initial condition where the dipole operator for the 80—7-82 tran- 
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sitions acts on monomer m = 1 only. Fig. |7] shows the subsequent population dynamics for two 
different values of the LVC strength. Focussing on the dimer (panel (a)) one notices that for small 
X exciton transfer between the two S2 excitation state dominates, while the Si state becomes pop- 
ulated only slowly. Upon increasing the LVC strength the local S2 state is rapidly deactivated and 
the transfer proceeds in the Si manifold. This behavior carries over to the trimer in Fig. Et,d, 
where the delayed population of states corresponding to excitations along the aggregate is even 
more visible. 

The most surprising result is perhaps the behavior of the Si state populations in the weak 
coupling case. Their population rise resembles that of the cases shown in Fig. |5t and e. This 
implies that the Si band is gradually filled by the rapid energy transfer in the S2 band without 
showing any signature of Si exciton transfer along the aggregate by it own. 

Summary 

Combining models from Frenkel exciton transfer and linear vibronic coupling theory the interplay 
of local deactivation processes and intermolecular energy transfer was investigated for a monomer 
as well as a dimer and trimer PBI system in dependence on the LVC strength. Although there is 
no obvious influence of LVC strength on low-resolution linear absorption spectra various features 
appeared in high-resolution spectra as well as in population and transfer dynamics, which can be 
explained by changes in the PES landscape due to the different couplings. Most notable here are 
the two opposite effects of the Coulomb coupling for a given LVC strength. First, the CIs are 
energetically shifted what may favor resonance with the optically excited wave packet. Second, in 
the Coulomb coupled case the wave packet does not leave the Franck-Condon region in direction 
of the CI, but only after backscattering from repulsive parts of the PES. The accompanying wave 
packet spreading diminishes the efficiency for non-adiabatic transitions at least in the weak LVC 
case. As net effect internal conversion becomes faster when going from the monomer to the dimer, 
but slows down for the trimer. It is interesting to note that a speed-up by a factors of about ten 



18 



upon dimerization has been observed in Ref. 
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although for a different system. 
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